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ABSTRACT 

We present the results of radio, optical and near-infrared observations of the field of 
TOC J0233. 3+3021, a cluster of milliJansky radiosources from the TexOx Cluster survey. In 
an observation of this field with the Ryle Telescope (RT) at 15 GHz, we measure a decrement 
in the cosmic microwave background (CMB) of — 675±95/iJy on the RT's w 0.65 kA baseline. 
Using optical and infrared imaging with the McDonald 2.7-m Smith Reflector, Calar Alto 3.5- 
m telescope and UKIRT, we identify the host galaxies of five of the radiosources and measure 
magnitudes of R w 24, J w 20, K w 18. 

The CMB decrement is consistent with the Sunyaev-Zel'dovich (SZ) effect of a massive 
cluster of galaxies, which if modelled as a spherical King profile of core radius 9q — 20" has a 
central temperature decrement of 900 //K. The magnitudes and colours of the galaxies are con- 
sistent with those of old ellipticals at z ~ 1. We therefore conclude that TOC J0233. 3+3021 
is a massive, high redshift cluster. These observations add to the growing evidence for a sig- 
nificant population of massive clusters at high redshift, and demonstrate the effectiveness of 
combining searches for AGN 'signposts' to clusters with the redshift-independence of the SZ 
effect. 
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1 INTRODUCTION 

Given a set of cosmological parameters such as Qm and Qa, N— 
body simulations of the build-up of the mass function of dark mat- 
ter halos (see e.g. Jenkins et al, 2000) are now sufficiently advanced 
that there are clear predictions concerning the number of systems 
of a given mass expected at a given cosmic epoch. At redshifts 
z « 1 and beyond, however, only a few members of the popula- 
tion of galaxy clusters are known (e.g. Stanford et al. 1997; Rosati 
et al. 1999; Chapman et al. 2000; Blanton et al. 2000; Fabian et al. 
2001). Here we present the first result of a new technique designed 
to find distant clusters. 

The TexOx Cluster (TOC) survey is an initiative designed 
to find massive high-redshift clusters. The TOC survey searches 
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for statistical overdensities of radiosources in the deep, wide area 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998). Although 
any search for clusters using AGNs as signposts will find only a 
small fraction of the total population, the surface density of ra- 
diosources is sufficiently low, and the sky coverage of the NVSS 
sufficiently large, that this technique is very efficient. TOC selects 
candidates where the radiosource overdensity is greater than 5-<r 
in a 7 x 7 arcmin 2 box. High-resolution VLA followup is used to 
eliminate the small number of cases where the apparent clustering 
of sources is spurious, and then optical and infrared imaging are 
carried out to identify the hosts galaxies of the radiosources and 
other candidate clusters galaxies. 

A full description of the first TOC results will be presented in 
a forthcoming paper (Croft et al. in preparation; hereafter C2001). 
In this letter we describe observations of one of the TOC candi- 
dates, TOC J0233. 3+3021. Our particular motivation in these ob- 
servations was to exploit the Sunyaev-Zel'dovich (SZ, Sunyaev and 
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Zel'dovich 1972) effect as a means of confirming candidate high- 
redshift clusters. 

In any cluster of galaxies, electrons in the intracluster plasma 
scatter the CMB photons passing through, distorting the thermal 
spectrum of the incident CMB and producing a change in bright- 
ness temperature given by ATsz oc J n e T e dl. Here n e and T e are 
the electron number density and Temperature, and dl is an element 
of line-of-sight element through the cluster. Two consequences of 
this are: 

• Because the CMB imprinting process is a scattering, the same 
ATsz is produced by a given J n e T B dl whatever the redshift of 
the cluster. 

• The SZ signal from a cluster is less sensitive to mass concen- 
tration than either X-ray (emissivity oc J iij dl) or optical searches. 
The integral of ATsz over the projected surface of a cluster is pro- 
portional to just the mass-temperature product of the gas. 

Observations of the SZ effect in several X-ray selected clus- 
ters at 2 > 0.5 have now been made using cm-wave interferometers 
(e.g. Grego et. al 2001, Joy et al. 2001, Grainge et al. 2001). Indeed, 
the SZ effect is an ideal tool for conducting a survey for clusters, 
not only because it allows cluster selection to high redshift, but 
because the selection is limited only by gas mass, almost indepen- 
dently of redshift. Plans have been developed for a new generation 
of dedicated SZ survey telescopes (e.g. Kneissl et al. 2001, Holder 
et al. 2000); however, in the absence of an SZ survey, we have un- 
dertaken pointed SZ observations of a TOC candidate to confirm 
that it is a high-z cluster and estimate its mass. 
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2 OBSERVATIONS AND RESULTS 
2.1 Ryle Telescope Observations 

The standard configuration of the Ryle Telescope (RT) for SZ ob- 
servations has described in detail elsewhere (e.g. Grainge et al. 
1996). A nearly-east- west array of five 13-m antennas is used, op- 
erating at 15 GHz, giving projected baselines from 13m to 108m 
(0.65-5.4 kA). Unresolved sources in the telescope field of view — 
typically radiogalaxies — contribute equal flux to all baselines of 
the array, whereas the extended structure of any SZ signal is typ- 
ically present only on the shortest baseline. We therefore use the 
long baselines to identify and remove point sources from the field, 
and use the residual flux on the shortest baseline to measure the SZ 
flux decrement. 

RT observations of the TOC J0233.3+3021 field totalling 280 
hours were made during 35 days in the period 1999 January to 1999 
June (the declination of the source means that it can only be ob- 
served for 8 hours each day). The telescope was in Cb configura- 
tion (Grainge et al. 1996), resulting in the aperture-plane coverage 
shown in Fig. []]. In this configuration the array provides two base- 
lines of maximum projected length 900A, three of length 1800A, 
two of 2700A and one each of 3600A, 4500A and 5400A. Note 
that since each day's observation is only eight hours, the baseline 
tracks in the uv plane are azimuthally truncated. For each day, ob- 
servations of the target field were interleaved with observations of 
a phase calibrator about every 20 minutes, and a primary flux cal- 
ibrator (3C286 or 3C48) was observed at either the start or end of 
the run. 

We analysed the combined 35-day visibility data for evidence 
of a flux decrement on the shortest baseline. First, we made and 
CLEANed a naturally-weighted map, using only visibilities from 



Figure 1. Aperture-plane coverage for the RT observations of the 
TOC J0233. 3+3021 field. Two physical baselines contribute to the shortest 
baseline in this figure, three to the second-shortest, two to the third shortest, 
and one to each of the three longest baselines. 



projected baselines longer than 1.5 kA, which is shown in Fig^, 
with the NVSS image of the TOC J0233.3+3021 field overlaid. 
Several of the NVSS sources are readily identifiable on the RT 
long-baseline map. Next, two of us (HJB and GC) independently 
investigated the sources present in the RT long-baseline map. We 
measured the positions of sources in the map using the task MAX- 
FIT in AlPS, then subtracted the brightest source from the visibil- 
ity data using UVSUB in AlPS. This process was repeated until no 
sources remained whose peak fluxes were greater than 3.5 times 
the rms in the map well outside the primary beam. A list of eight 
sources was agreed upon; all but one of these sources is identifi- 
able on the NVSS map, as shown in Fig|^. Accurate positions and 
fluxes for these sources were then calculated using the the FLUX- 
FlTTER algorithm (Grainger et al. 2001). FLUXFlTTER takes input 
estimates of the point-source fluxes and positions, and varies these 
to minimize x-squared between the visibility data and model visi- 
bilities containing the point sources and noise; it has been demon- 
strated to produce robust flux and position measurements (Grainger 
et al. 2001). The resulting source positions and fluxes, calculated 
using all the data from the projected baselines longer than 1.5 kA 
shown in Table 2. 1 . These sources were then subtracted from all the 



visibilites. To check that no residual flux from any of the sources 
was present, source-subtracted maps were made using baselines of 
length 1.5 kA and greater (shown in Fig ^|), and 2.0 kA and greater. 
These maps are both consistent with noise. 

Finally, a map was made using source-subtracted visibilities 
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Table 1. Details of the radiosources in the TOC J0233. 3+3021 field, in the following format: Position of source at 15-GHz as subtracted from RT observations; 
flux at 15 GHz; Position of host galaxy identification (C2001) ; R-, J- and i^-band magnitudes (C2001); Label used in Fig. 2. 



RTRA(J2000) RTDec(J2000) Si 5G Hz ID RA (J2000) ID Dec (J2000) R J K Label 

/mJy 

0233 18.22 3022 42.91 2.007 0233 18.1 302240.1 23.75 19.75 18.19 A 

02 33 24.97 302009.31 1.906 02 33 25.0 3020 11.2 23.28 19.74 > 19 B 

02 33 11.41 3025 07.91 1.118 

02 33 31.17 30 22 33.31 1.077 02 33 31.3 30 22 31.2 23.35 20.76 18.58 C 

02 33 23.65 30 22 06.41 1.053 02 33 23.7 30 22 00.9 23.10 19.86 18.19 D 

02 33 14.23 30 19 46.41 0.497 02 33 14.1 30 19 39.8 24.96 20.82 18.39 E 

02 33 38.87 30 21 10.41 0.352 

02 33 22.83 30 19 24.11 0.339 
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Figure 2. CLEANed 1.5-5.4kA Ryle 

Telescope map of the TOC J0233. 3+3021 field before source subtraction 
overlaid with the NVSS 1.4-GHz map of the field (contours). The greyscale 
runs from —0.5 to 0.6 mJy beam -1 , and the contours are 3, 3\/2, 6, 6\/2, 
. . . multiples of 0.6 mJy beam -1 . The large circle shows the FWHM of the 
RT primary beam, and the ellipse at lower left shows the FWHM of the RT 
synthesized beam. The five NVSS sources which define the TOC survey 
overdensity are labelled A — E; the crosses show the positions of th e po int 
sources found at 15 GHz. Details of these sources are listed in Table |2.l|. 



with projected baselines in the range 0.65-1.0 kA, corresponding 
to the single shortest physical baseline. This map showed a nega- 
tive feature centred at 02 h 33 m 20.8 s , 30° 21' 58.7" (J2000), and 
with a minimum flux of in the dirty map of — 675/iJy. The syn- 
thesized beam of this single-baseline map has strong sidelobes, so 
the rms level on the map even well outside the primary beam is 
still somewhat affected by the presence of the negative feature at 
the centre of the map. We therefore made a very large map (2048 
x 2048 5"-square pixels) and measured the map rms in 128-pixel 
square regions where the map apppeared by eye to be unaffected 
by sidelobes from the central source. The rms levels agreed wth 
the expected noise, due to the telescope system temperature, of 
95 /iJy. The short-baseline map was CLEANED to a flux limit of 
— 285/iJy, and this map is presented in Fig. El On inspection of the 
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Figure 3. 1.5-5.4kA Ryle Telescope map of the TOC J0233.3+3021 field 
after subtraction of the point sources listed in Table |2.l| . The greyscale runs 
from —0.217 to 0.296 mJy beam -1 , and contours are plotted at -3 and 3 
times 0.050 mJy beam -1 , the map rms outside the primary beam. The map 
is consistent with noise. The primary and synthesized beams are shown as 
in Fig. ^. The small circles are positions of the five radiogalaxy optical IDs 
given in Table 1 . 



azimuthally-averaged visibilities (Fig. |5|), the flux on the shortest 
baseline is — 695/^Jy, with an rms of 170 /iJy (note that since the 
shortest baseline has an elliptical track in the aperture plane, this 
rms includes the variation in the amplitude of the SZ signal as a 
function of uv distance). 

2.2 Optical and infra-red observations of the galaxies 

We obtained .R-band imaging of the TOC J0233.3+3021 field us- 
ing the Imaging Grism Instrument on the McDonald Observatory 
2.7-m Smith Reflector in 1998 October, 1999 October, and 2000 
September, and J-band imaging of the field at the Calar Alto 3.5-m 
telescope in 1999 August. Additionally, we obtained several K- 
band images of the TOC J0233. 3+3021 field with UFTI at UKIRT 
in 2000 July. Full details, and extensive discussion of other candi- 
date cluster galaxies, will be presented in C2001. Using these im- 



© 0000 RAS, MNRAS 000, 000-000 



4 



Garret Cotter et al. 




RIGHT ASCENSION (J2000) 



Figure 4. CLEANed 0.65-1.0 

kA Ryle Telescope map of the TO C J0 233. 3+3021 field after subtraction 
of the point sources listed in Table |2~l| . The greyscale runs from —0.675 
to 0.335 mJy beam - 1 . Contours are plotted at -7,-6,-5,-3,-4,-2 and 2 times 
0.095 mJy beam -1 , whihc is consistent with the map rms well outside the 
primary beam, as described in Section 3.2. Beams and source positions are 
shown as in previous figures. 
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Figure 5. Real (top) and imaginary (bottom) parts of the visibilities, az- 
imuthally averaged about the centre of the decrement after source subtrac- 
tion. The baseline scale ranges from to 6000 wavelengths, and the flux 
scale in each plot ranges from —1.0 to 1.0 mJy. The visibilities indicate a 
single extended negative source at the centre of the field. 



ages, in combination with high-resolution VLA imaging (C2001), 
we have identified the host galaxies of the central five NVSS ra- 
diosources shown in Fig. ^. For four of these galaxies, we have 
obtained 3.5"-diameter aperture magnitudes in all filters; for the 
fifth galaxy, only R and J magnitudes are available. The measured 
magnitudes of the galaxies are presented in Table 2.1 along with 



the positions of the host galaxy identifications. The positions were 
measured from the i?-band image, with an astrometric solution us- 
ing the positions of stars measured in the Cambridge APM survey. 
The magnitudes of the galaxies are all roughly K w 18, J « 20, 
R » 24. 

In 2000 January, we attempted spectroscopy of the radiogalax- 
ies with ISIS on the WHT. We used a 2-arcsec slit, with sensitivity 
on the ISIS red arm up to 9340 A. We took 1800-second spectra of 
A (at PA=1 18, the axis of its double structure; see C2001), galaxies 
B & C (as a pair with PA=1 13), and galaxies D & E (as a pair with 
PA=77). Faint red continuum was detected spectroscopically in all 
the cases expected (i.e. the 4 brightest), but no emission lines were 
evident. 



3 DISCUSSION 

3.1 The redshifts of the radiogalaxies 

The optical and infrared magnitudes of the radiogalaxies immedi- 
ately imply that they lie at significant redshift. In particular, given 
the if-band Hubble diagram for faint radiosources (such as the 7C 
sample of Willott, Rawlings and Blundell, 2001), the K ~ 18 mag- 
nitudes of the TOC J0233. 3+3021 radiogalaxies suggest that they 
are L* galaxies at redshifts of z ~ 1.0. We next consider the RJK 
colours of the radiogalaxies. It is now well-established that faint ra- 
diogalaxies at z ~ 1 are hosted by Extremely Red Objects (EROs) 
with colours as red as R — K » 6. (e.g. Dunlop et al. 1996; Willott, 
Rawlings and Blundell, 2001). This is consistent with a picure in 
which faint radiosources at these redshifts are hosted by old ellip- 
tical galaxies. Comparing the the radiogalaxies with galaxies in the 
z = 1.27 cluster C1G J0848+4453 found by Stanford et al. (1997), 
we find a close similarity in RJK colour-magnitude diagrams (Fig. 
^|). If we plot the RJK colours of the radiogalaxies against the 
track of a redshifted E0 galaxy, we find that they cluster about the 
z » 1 region on colour-colour space (Fig. Q). Lastly, the faint red 
continuum detected in the optical spectra of the radiogalaxies is 
also consistent with their lying at z w 1. 

These optical and infrared data do not alone demonstrate un- 
equivocally that this overdensity of radiogalaxies is a real cluster. 
However, the coincidence between the overdensity of sources and 
the CMB decrement leads us to conclude that we are seeing the SZ 
effect of a massive cluster, and that most, if not all, the radiogalax- 
ies are cluster members. 



3.2 Cluster properties inferred from the SZ effect 

From the SZ flux on the shortest spacing, we can now estimate the 
gas mass of the cluster. We used PROFILE (Grainge et al. 2001b) to 
model the cluster as a spherical King-profile gas distribution, with 
a gas temperature of 7keV, /3-parameter 0.65, and find the best fit 
value of the central electron density (no) to the data for varying core 
radius 9c- While the temperature and /3 parameter are unknown 
in the absence of deep X-ray imaging, the values we assume are 
typical of those which have been been measured in X-ray selected 
z > 0.5 clusters (see, e.g., Grego et al. 2001 for a summary). For 
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Figure 6. R — K and J — K colour-magnitude plots for the host galaxies 
of the radiosources; cf Fig. 4 of Stanford et al. (1997). The dashed lines 
show the red sequence of the galaxies in the z = 1.27 cluster identified by 
Stanford et al. (1997). The similar red colours of the TOC J0233.3+3021 
galaxies suggests that they are old elliptical galaxies, with the R and J 
filters straddling the observed-frame 4000-A break. 
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an Q.m = 0.3, Qa = 0.7 cosmology, with Ho = 65kms 1 Mpc 1 , 
we find that the gas mass within a radius of 1 Mpc is proportional 
to Oq. We place the constraint 6c > 20" from our measurement 
that the flux on the second shortest RT baseline is less than 200 ujy. 
This core radius corresponds to a physical distance of 178 kpc and 
the appropriate best fit model has no = 10 4 m~ 3 and a central 
decrement of 900 fiK. We therefore estimate that the minimum gas 
mass of the cluster is 5 x 10 13 Mq. Assuming that the total mass 
of the cluster is ten times greater than this, TOC J0233. 3+3021 is 
clearly a massive cluster similar to those identified at z > 0.5 in 
X-ray surveys such as EMSS. 

The accurate measurement of the properties of this population 
of high- 2 massive clusters is a significant goal in observational cos- 
mology. In this Letter, we have demonstrated that current SZ obser- 
vations are effective in confirming high-z cluster candidates which 
have been selected by other means; SZ investigations of other TOC 
candidates are continuing. The new generation of SZ survey tele- 
scopes will not be restricted to the use of 'signposts' to high-2 
clusters, and promise to detect the majority of the distant cluster 
population. 



4 CONCLUSIONS 

We have observed the field of TOC J0233.3+3021, a cluster of mil- 
liJansky radiosources selected from the NRAO VLA Sky Survey 
(NVSS), with the Ryle Telescope (RT) at 15 GHz, and in the opti- 
cal and near-infrared with the McDonald 2.7-m, Calar Alto 3.5-m, 
and UKIRT. Our conclusions are as follows: 

(i) After subtraction of the baseline-independent fluxes of the 
radiosources, we find a seven-sigma decrement on the RT shortest 
baseline of —695 /xJy. 

(ii) The radiosource host galaxies have magnitudes of K « 18, 
J ~ 20, R S3 24. From the infrared magnitude-redshift relation for 
radiosources, the galaxies must lie at z > 0.5. The similar RJK 
colours strongly suggest that these galaxies are evolved ellipticals 
at z ~, and the galaxy magnitudes are consistent with the R — J 
colour-magnitude diagrams of known z > 1 clusters. 

(iii) We therefore conclude that we have identified a strong SZ 
decrement with a cluster of galaxies at z ~ 1. To produce the SZ 
decrement, the cluster must have a mass-temperature product sim- 
ilar to that of massive Abell clusters at low redshift; we estimate a 
minumum gas mass of 5 x 10 13 M Q . These observations support 
the growing evidence for a population of massive clusters at z > 1. 

(iv) These observations demonstrate the effectiveness of com- 
bining searches for AGN 'signposts' to clusters with the redshift- 
independence of the SZ effect. Investigations of this and other high- 
z cluster candidates selected from NVSS are continuing. 



Figure 7. The host galaxies of the radiosources in the RJK colour-colour 
plane. The radiogalaxies are shown as large crosses, with the size of the 
cross indicting the photometric erros. The line shows the colours of a non- 
evolving E0 galaxy as it is redshifted from z = 0.0 (lower left), with sym- 
bols marking intervals of Az = 0.1. The centroid of the radiosources in 
colour-colour space is at z S3 1. 
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